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We present recent results on lattice simulations using chiral effective field theory. In particular 
we discuss lattice simulations for dilute neutron matter at next-to-leading order and three-body 
forces in light nuclei at next-to-next-to-leading order. 
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Lattice simulations based on the framework of effective field theory have been used in studies 
of nuclear matter [|l|] and neutron matter [|2|, ^, The method has also been used to study light 
nuclei in pionless effective field theory [||] and chiral effective field theory at leading order (LO) 

More recently next-to-leading order (NLO) calculations have been carried out for the ground 
state of neutron matter A review of lattice effective field theory calculations can be found 

in Ref. In this proceedings article we describe some new results which, at the time of this 
writing, are not yet published. 

At leading order in chiral effective field theory the nucleon-nucleon effective potential is 

VlO=V+V;2+VOPE^ (1) 

V , Vj2 are the two independent contact interactions at leading order in the Weinberg power counting 
scheme, and \/0PEP ^j^g instantaneous one-pion exchange potential. The interactions in Vlo can 
be described in terms of their matrix elements with two-nucleon incoming and outgoing momentum 
states. For bookkeeping purposes we label the amplitude as though the two interacting nucleons 
were distinguishable, A and B. In the following q denotes the ^-channel momentum transfer while 
k is the M-channel exchanged momentum transfer. We use T/ with / = 1,2,3 to represent Pauli 
matrices acting in isospin space and Os with S = 1,2,3 to represent Pauli matrices acting in spin 
space. 

For the two leading-order contact interactions the amplitudes are 

£/{V)=C, (2) 
^iVp)=Cj2l^TfTf. (3) 

For the one-pion exchange potential. 



For our physical constants we take m = 938.92 MeV as the nucleon mass, mji = 138.08 MeV as 
the pion mass, fji = 93 MeV as the pion decay constant, and gA = 1.26 as the nucleon axial charge. 

In Ref. two different lattice actions were considered which were later denoted LOi and LO2 
[^]. The interactions in VlOi include one-pion exchange and two zero-range contact interactions 
corresponding with amplitude 

^ ( Vlo, )=C + Ct2Y.xtzf + s^ (y OPEP) . (5) 

The interactions in VLO2 consist of one-pion exchange and two Gaussian-smeared contact interac- 
tions, 

^ ( Vlo, ) = Cf{q) + C,2m £ rf xf + ^ (yOPEP) , (6) 

/ 

where f{q) is a lattice approximation to a Gaussian function. The smeared interactions in LO2 
are used to better reproduce ^-wave phase shifts for nucleon momenta up to the pion mass. The 
coefficients C and C/2 are tuned to reproduce the physical S'-wave scattering lengths. In Ref. [^] 
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nucleon-nucleon phase shifts were calculated for these two lattice actions using the spherical wall 



method [10| at spatial lattice spacing a = (100 MeV)^^ and temporal lattice spacing at = (70 
MeV)^^. For each case NLO corrections were also computed perturbatively and the unknown 
operator coefficients determined by fitting to low-energy scattering data. 

In Ref. the ground state energy for dilute neutron matter was computed using the lattice 
action LO2 and auxiliary-field Monte Carlo. Next-to-leading-order corrections to the energy were 
also calculated perturbatively. In this calculation the largest source of systematic error was the 
large size of NLO corrections for Fermi momenta larger than 100 MeV. This was due to attractive 
P-wave interactions generated by Gaussian smearing in LO2 that needed to be cancelled at next- 
to-leading order. In systems with both protons and neutrons this P-wave correction is numerically 
small when compared with the strong binding produced by 5-wave interactions. For pure neu- 
tron matter, however, the 5-wave interactions produce much less binding due to Fermi repulsion. 
Therefore on a relative scale, the P-wave interactions are not as small an effect in neutron matter. 

These problems have been resolved using a new leading-order action LO3 [|^]. The interac- 
tions in Vlos correspond with the amplitude, 

'1 1^ . «\ /3 1 

The Gaussian-smeared interactions are multiplied by spin and isospin projection operators. Only 
the Cs=o,/=i term contributes in pure neutron matter. Using the LO3 action with NLO conections, 
we have computed the ground state energy for dilute neutrons in a periodic box [^]. For spatial 
lattice spacing a = (100 MeV)^' and temporal lattice spacing at = (70 MeV)^' simulations were 
done with 8, 12, 16 neutrons in periodic boxes with lengths L = 4,5,6,7. In Fig. [l]we show results 
for the ratio of the interacting ground state energy to non-interacting ground state energy, Eo/E^^^, 
as a function of Fermi momentum kp. For comparison we show other results from the literature: 
FP 1981 [|ri|], APR 1998 [|T2|], CMPR v6 and v8' Q, SP 2005 GC 2007 and GIFPS 



^ (Vlo, ) = Cs=OJ=J{q) 4 - 4 I C^l 4 + 4 I ^f^f 



2008 ||T§]. 

At next-to-next-to-leading order (NNLO) in chiral effective field theory we find contributions 
due to three-nucleon forces. These interactions consist of a pure contact interaction, one-pion 
exchange, and two-pion exchange [p^]. The coupling of one or more pions to a single nucleon is 



constrained by chiral symmetry and the corresponding low energy constants are known [|18|]. In the 
limit of exact isospin symmetry there are only two unknown coefficients, one for the three-nucleon 
contact interaction and one for the two-nucleon-pion vertex involved in the one-pion exchange 
interaction. At fixed lattice spacing we have determined these two unknown coefficients by fitting 
to the triton binding energy and spin-doublet nucleon-deuteron scattering phase shifts via Liischer's 
finite volume formula [[1^]. Results for the doublet nucleon-deuteron scattering phase shift are 
shown in Fig. ^ using the LO2 lattice action for lattice spacing a = (100 MeV)^' and temporal 
lattice spacing a, = (150 MeV)"^ [^]. 

Having determined the NNLO three-body forces, we have computed the ground state of the 
alpha particle without Coulomb interactions on a periodic lattice using auxiliary-field projection 
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Figure 1: Ground state energy ratio Eq/E^^^ for LO3 and NLO3 versus Fermi momentum kf- For com- 
parison we show results for FP 1981 APR 1998 CMPR v6 and v8' 2003 SP 2005 [|l4[, GC 
2007 and GIFPS 2008 [|l6|. 
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Figure 2: Results for the spin-doublet nucleon-deuteron scattering phase shift at LO, NLO, and NNLO. 



Monte Carlo [|20|]. The NNLO results are within 5% of the actual Coulomb-subtracted alpha bind- 
ing energy of about 29 MeV. This is consistent with the expected size of errors for our chosen 
lattice spacing and order in effective field theory. 
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